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.  Research  Objectives 


The  research  work  on  this  project  has  been  concerned  with 
dynamic  effects  on  fracture.  Two  main  areas  have  been  investigated: 
"high-rate  loads  on  bodies  containing  cracks."  and  "fast  fracture 
and  crack  arrest." 

Dynamic  effects  become  important  if  the  external  loads  give  rise 
to  propagating  mechanical  disturbances  (as  for  impact  loads  and 
explosive  charges)  which  can  strike  a  crack  and  cause  crack  propagation. 
Spalling  is  an  example  of  a  fracture  phenomenon  caused  by  the  rapid 
application  of  loads.  Dynamic  effects  become  also  important  if  a 
crack  propagates  very  rapidly,  so  that  rapid  motions  are  generated  in 
the  solid.  Even  though  in  most  cases  fracture  studies  should  focus  on 
the  prevention  of  crack  propagation,  cracks  sometimes  become  unstable 
despite  the  best  attempts  for  prevention.  Both  for  brittle  and 
ductile  materials,  it  is  then  important  to  know  whether  high  speeds  of 
crack  propagation  can  be  expected,  what  the  dynamic  Interaction  with 
the  geometry  of  the  body  will  be,  and  most  Importantly  how  a  rapidly 
propagating  crack  can  be  arrested. 

A  summary  of  the  research  objectivs  of  the  completed  project  is 
DYNAMIC  EFFECTS  ON  FRACTURE 

I.  HIGH-RATE  LOADS  ON  BODIES  II.  FAST  FRACTURE  AND  CRACK  ARREST. 

CONTAINING  CRACKS. 


I-A:  Computation  of  elastodynamic 
stress  intensity  factors. 

I-B:  Influence  of  plastic  yielding 
on  the  dynamic  stress  field 
near  a  crack  tip. 

I-C:  Conditions  for  crack  propaga¬ 
tion  under  dynamic  loading 
conditions. 


II-A: Computation  of  elastodynamic 
stress  intensity  factors  for 
rapidly  propagating  cracks. 

II-B: Analysis  of  plastic  deforma¬ 
tions  near  a  rapidly  propagating 
crack  tip. 

II-C:Application  of  fracture 

criteria  to  compute  crack-tip 
speeds . 

II-D:Arrest  of  rapidly  propagating 
cracks . 


II -E: Crack  forking  and  kinking 

of  rapidly  propagating  cracks 


The  work  has  evolved  through  four  stages  which  were  stated 
in  the  original  proposal.  These  stages  are 

Stage  1:  Development  of  analytical  and  numerical  methods  to 
analyze  dynamic  effects  on  crack-tip  fields. 

Stage  2:  Comparison  of  crack-tip  fields  for  linear  elasticity 
and  various  constitutive  models  of  elasto-plastic  and  visco-plastic 
material  behavior. 

Stage  3:  Comparison  of  various  fracture  criteria. 

Stage  4;  Application  of  fracture  criteria  and  the  formulation  of 
conditions  for  initiation  of  crack  propagation,  continued  crack 
propagation  and  crack  arrest,  as  well  as  for  crack  kinking  and  crack 
forking  under  dynamic  conditions. 

From  the  point  of  view  of  the  Air  Force,  the  technological  signi¬ 
ficance  and  relevance  of  the  completed  research  derives  from  the  fact 
that  it  is  at  the  present  time  either  technically  or  economically  not 
feasible  to  manufacture  so  that  cracks  never  occur  in  factory-new 
structures,  or  after  a  period  of  service.  It  is,  therefore,  necessary 
to  be  able  to  estimate  the  stability  of  flaws,  the  propagating 
characteristics  for  fast  fracture  and  the  arrest  mechanisms  for  rapidly 
propagating  cracks,  under  a  variety  of  loading  conditions  which  should 
include  conditions  generated  by  Impact  loads  and  explosive  charges. 


2.  Summary  of  Completed  Research 


Journals  and  Proceedings  of  Conferences.  A  complete  list  of  publications 
is  given  in  Appendix  A.  Lectures  and  presentations  of  material  completed 
on  this  Project  are  listed  in  Appendix  B. 


The  research  results  are  briefly  summarized  as  follows. 

2.1  Effects  of  Plasticity  on  Near-Tip  Fields 

Under  rapid  loading  conditions  the  mass  density  of  a  material 

affects  the  fields  of  stress  and  deformation  near  a  crack  tip.  For  such 

dynamic  problems  plastic  deformation  in  the  immediate  vicinity  of  a  crack 

* 

tip  was  investigated  in  [1]  for  a  stationary  crack.  Deformation  theory 
was  employed  for  the  first  phase  of  the  loading  when  the  fields  are 
increasing  monotonically  with  time.  The  general  character  of  the  near-tip 
fields  was  analyzed  both  with  respect  to  its  variation  with  time  and  with 
polar  angle.  The  non-linear  near-tip  fields  were  related  to  the  linearly 
elastic  far-field  by  means  of  a  path-independent  Integral. 

In  [2],  deformation  theory  was  applied  to  a  thin  plate  containing  a 
through -crack.  The  in-plane  normal  stresses  were  assumed  to  be  tensile, 
but  with  a  non-uniform  variation  through  the  thickness  of  the  plate,  as 
a  result  of  bending  and  a  dominant  extension.  The  nonlinear  near-tip 
fields  (which  are  singular)  have  been  analyzed  asymptotically  on  the  basis 
of  plate  theory.  It  was  found  that  the  angular  variations  of  the  near-tip 
fields  are  just  the  same  as  for  generalized  plane  stress.  Assuming  small- 
scale  yielding  a  path-independent  Integral,  which  is  valid  in  a  region 
close  to  the  crack  edge,  was  used  to  connect  the  nonlinear  near-tip  fields 
with  the  corresponding  singular  parts  of  the  linear  fields.  It  was  shown 
that  thertmlinear  behavior  near  the  crack  edge  significantly  affects  the 
through-the-thickness  variations  of  the  near-tip  fields.  The  singular 
parts  of  the  in-plane  stresses  tend  to  become  more  uniform  through  the 
thickness  of  the  plate  with  stronger  strain  hardening. 


* 


Numbers  in  brackets  refer  to  the  list  of  publications  in  Appendix  A. 


For  dynamic  fracture  problems  plastic  deformation  in  the  immediate 
vicinity  of  a  crack-tip  was  investigated  in  [3].  Both  stationary  and 
propagating  crack  tips  were  considered.  For  a  stationary  crack  tip, 
deformation  theory  was  employed  for  the  first  pha^e  of  the  loading  when 
the  fields  are  increasing  monotonically  with  time.  The  general  character 
of  the  near-tip  fields  was  analyzed  both  with  respect  to  its  variation 
with  time  and  with  polar  angle.  The  non-linear  near-tip  fields  were 
related  to  the  linearly  elastic  far-field  by  means  of  a  path-independent 
integral.  In  the  second  part  of  the  paper  we  considered  rapidly  propa¬ 
gating  cracks.  The  near-tip  fields  for  various  models  of  material 
behavior  were  discussed.  In  particular,  some  earlier  work  by  Achenbach 
and  Kanninen  for  a  rapidly  propagating  Mode-Ill  crack,  in  a  material  which 
displays  strain  hardening  was  reviewed.  In  the  last  part  of  the  paper  the 
fields  near  a  rapidly  propagating  crack-tip  in  an  elastic  perfectly-plastic 
material  was  considered  for  the  case  that  inertia  terms  are  of  importance. 
The  system  of  governing  equations  in  the  plastic  region  was  presented  and 
shown  to  be  hyperbolic  in  nature.  As  a  first  approximation  the  steady- 
state  case  with  respect  to  the  moving  crack-tip  was  considered  and  an 
asymptotic  analysis  of  the  near-tip  field  was  carried  out. 

A  dynamic  analysis  of  fast  fracture  and  crack  arrest  based  on  the 
Dugdale  model  was  presented  in  [7].  Numerical  approaches  by  the  finite 
difference  method  have  been  considered  in  [13]  and  [18] . 

2.2  Computation  of  Elastodynamic  Stress-Intensity  Factors 

The  effect  of  proximity  of  a  boundary  on  elastodynamic  stress 
intensity  factors  has  been  investigated  in  [4]  for  a  surface-breaking 
crack,  and  in  [14]  for  a  subsurface  crack  parallel  to  the  surface  of  a 
half-space.  In  recent  work  we  have  also  considered  a  subsurface  crack 
which  is  oriented  under  an  arbitrary  angle  with  the  surface  of  a  half¬ 
space,  see  [16]. 


The  configurations  that  have  been  considered,  are  two  dimensional 
with  deformations  in  plane  strain.  Systems  of  coupled  singular  integral 
equations  for  the  Mode-I  and  Mode-II  dislocation  densities  have  been 
derived.  These  equations  have  been  solved  numerically  for  the  cases  of 
time-harmonic  uniform  tension  and  uniform  shear  applied  at  the  surface 
of  the  half-space.  The  ratios  of  elastodynamic  to  eleastostatic  stress 
intensity  factors  have  been  computed.  The  results  display  the  dependence 
on  the  frequency  and  on  the  ratio  d/a,  where  d  is  the  distance  from  the 
upper  crack  tip  to  the  free  surface  and  a  is  the  crack  length.  For 
small  angles  of  inclination  with  the  free  surface,  and  for  small  values 
of  d/a,  time-harmonic  excitations  of  the  body  may  induce  quite  strong 
resonance  vibrations  of  the  layer  between  the  crack  and  the  free  surface. 
Such  resonance  vibrations  give  rise  to  substantial  increases  in  both  the 
Mode-I  and  Mode-II  stress  intensity  factors.  These  resonance  effects 
were  Investigated  in  some  detail  for  the  parallel  crack  in  [14]. 

A  three-dimensional  stress  analysis  problem  for  a  surface-breaking 
crack  was  considered  in  [17].  A  half-space  containing  a  surface-breaking 
crack  of  uniform  depth,  oriented  normal  to  the  surface,  was  subjected  to 
three-dimensional  dynamic  loading.  The  elastodynamic  stress  analysis 
problem  was  decomposed  into  two  problems  which  are  symmetric  and 
antisymmetric,  respectively,  relative  to  the  plane  of  the  crack.  The 
formulation  of  these  problems  was  subsequently  reduced  to  singular 
integral  equations  for  the  corresponding  dislocation  densities,  namely, 
a  single  integral  equation  for  the  symmetric  problem  and  a  set  of 
coupled  integral  equations  for  the  antisymmetric  problem.  The  systems 
of  integral  equations  were  solved  numerically.  The  dislocation  densities 
directly  yield  the  stress-intensity  factors.  As  an  example  an  applied 
stress  field  with  harmonic  variation  along  the  crack  length  and  exponen¬ 
tial  decay  with  crack  depth  was  considered.  The  dependence  of  the  stress 
Intensity  factors  on  the  frequency  and  on  the  wavelength  in  the  crack 
direction  has  been  displayed  by  numerical  results . 


In  this  work,  we  have  investigated  both  the  effects  of  plastic 
deformation  near  a  propagating  crack  tip  and  dynamic  effects  due  to  high 
crack-tip  speeds.  The  constitutive  equations  that  have  been  employed 
define  an  elastic  viscoplastic  material.  The  constitutive  model,  which 
was  proposed  by  Bodner  and  Partom  does  not  require  the  statement  of  a 
separate  yield  criterion,  nor  is  it  necessary  to  consider  loading  and 
unloading  separately.  Plastic  deformations  always  exist,  but  they  are 
negligibly  small  when  the  material  behavior  should  be  essentially  elastic. 

The  geometry  that  has  been  considered  is  a  two-dimensional  one  of  a 
thick  strip  which  contains  a  rapidly  propagating  semi-infinite  crack  in 
its  center-plane.  The  faces  of  the  strip  were  subjected  to  uniform  in¬ 
plane  displacements,  so  that  the  crack  propagates  in  Mode-I.  A  steady- 
state  situation  relative  to  the  moving  crack  tip  has  been  assumed.  The 
plastic  deformations  near  the  crack  tip,  the  residual  plastic  strains  in 
the  wake  of  the  crack  tip  and  other  field  variables  have  been  obtained 
directly  from  the  complete  solution  without  the  assumption  of  small  scale 
yielding. 

Detailed  results  have  been  given  in  [5] ,  [6]  and  [11].  The  method  of 
solution  is  numerical  and  involves  an  iterative  procedure  which  is 
continued  until  a  steady  state  solution  is  reached,  in  which  the 
equation  of  motion,  the  flow  rules  and  the  boundary  conditions  are 
satisfied  simultaneously.  The  method  is  based  on  a  finite  difference 
procedure  which  is  unconditionally  stable  and  has  a  second  order  accuracy. 

In  the  special  case  of  a  perfectly  elastic  strip  the  problem 
possesses  an  analytical  solution  which  has  been  employed  to  check  the 
accuracy  of  the  numerical  method,  and  excellent  agreement  has  been 
obtained. 

Typical  effects  due  to  viscoplastic  constitutive  behavior  have  been 
studied  by  comparisons  with  the  corresponding  elastic  fields.  The  effects 
of  high  crack-tip  speeds,  which  are  directly  related  to  the  strain-rate 


dependence  of  the  material,  have  been  studied  by  comparisons  of  solutions 
for  three  crack-tip  velocities.  The  influence  of  the  inertia  term  in  the 
governing  equations  has  been  studied  by  comparisons  with  the  corresponding 
quasi-static  solutions.  In  particular,  the  dependence  on  the  crack-tip 
speed  of  the  plastic  zone  in  the  vicinity  of  the  crack  tip,  the  level  of 
plastic  straining,  the  amount  of  dissipative  plastic  work  and  the  crack¬ 
opening  displacements  have  been  examined. 

In  [9]  and  [12]  the  transient  problem  of  deceleration  and  arrest  of 
a  rapidly  propagating  crack  tip  has  been  investigated  for  a  crack  which 
initially  propagates  in  an  elastic  solid  but  then  enters  a  region  of 
viscoplastic  material  properties,  of  particular  interest  in  the  dynamic 
arrest  problem  is  the  plastic  strain  just  ahead  of  the  arresting  crack 
tip.  This  strain  has  been  computed  as  a  function  of  time.  Within  the 
context  of  a  critical  strain  criterion  it  has  been  assumed  that  the 
crack  tip  will  arrest  if  this  strain  does  not  exceed  a  certain  critical 
value . 

2.4  Crack-Kinking  under  Stress-Wave  Loading 

When  a  stress  pulse  strikes  a  crack,  the  crack  may  be  induced  to 
propagate,  but  not  necessarily  in  it's  own  plane.  In  earlier  work, 
it  has  been  attempted  to  explain  kinking  of  a  crack  at  finite  kinking 
angles.  In  work  completed  on  this  project,  [10]  and  [15],  we  have 
reconsidered  the  two-dimensional  configuration  of  an  initially  stationary 
crack  which  kinks  under  an  angle  ktt  with  it's  original  plane.  We  have 
discussed  earlier  results  for  Mode  III,  and  we  have  proposed  a  way  to 
approximate  the  elastodynamic  stress  intensity  factors  for  the  Mixed 
Mode  I-II  case. 

The  tip  of  the  kinked  crack  is  assumed  to  propagate  at  a  constant 
velocity  c^,,  and  kinking  is  initiated  at  an  angle  kit  at  the  instant 
that  an  incident  disturbance  first  strikes  the  original  crack  tip. 

These  two  assumptions  render  the  solution  self-similar.  Three  cases  have 
been  considered  corresponding  to  incidence  of  either  an  anti-plane 
transverse  wave,  an  in-plane  transverse  wave  or  a  longitudinal  wave. 
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The  elastodynamic  stress  intensity  factors  have  been  computed  as 
functions  of  the  crack  tip  speed,  c  ,  the  kinking  angle,  ktt,  and  the 
angle  of  wave  incidence  air.  For  a  given  angle  of  incidence,  the 
elastodynamic  stress  intensity  factors  have  been  used  to  compute 
the  corresponding  energy  fluxes  into  the  propagating  crack  tip. 

Mode-Ill  problems  of  the  kind  formulated  in  [10]  can  be  solved 
rigorously.-  Results  for  Mode-Ill  kinking  under  an  arbitrary  angle 
were  given  earlier.  The  corresponding  mixed  Mode  I-II  problems  have, 
however,  as  yet  eluded  a  rigorous  analytical  solution. 

For  the  mixed  Mode  I-II  problems  we  have  proposed  a  simple 
approximation  to  the  solution  of  the  superposition  problem.  The 
approximation  is  based  on  an  observation  from  the  exact  Mode-Ill 
solution  of  the  analogous  superposition  problem,  that  for  an  important 
range  of  kinking  angles  the  elastodynamic  stress  intensity  factor  of 
the  kinked  crack  is  affected  more  by  the  dependence  on  k  of  the  loading 
on  the  new  crack  faces  than  by  the  wedge  geometry  at  the  original  crack 
tip.  This  observation  then  suggests  that  in  first  approximation  we  may 
ignore  the  wedge  geometry  altogether,  and  we  may  compute  elastodynamic 
stress  intensity  factors  by  considering  a  crack  propagating  in  its  own 
plane  but  where  the  new  crack  faces  are  subjected  to  tractions  corres¬ 
ponding  to  the  kinking  crack. 

For  the  Mode-Ill  case  the  approximation  of  the  elastodynamic 
stress  intensity  factor  for  the  kinked  crack  can  be  checked  by 
comparisons  with  exact  results.  The  range  of  kinking  angles  Kir  for 
which  the  approximation  gives  good  results  turns  out  to  be  sur¬ 
prisingly  large.  For  the  mixed  Mode  I-II  case  comparisons  with  the 
numerical  results  have  been  carried  out,  and  satisfactory  agreement 
has  been  obtained.  No  rigorous  mathematical  proof  of  the  approxi¬ 
mation's  validity  is  given  in  this  paper.  It  is,  however,  to  be 


expected  that  the  results  correspond  to  the  first  order  terms  in  a 
perturbation  procedure  for  small  kinking  angle. 

The  results  of  this  work  suggested  that  the  approximation  could 
be  used  to  analyze  crack  kinking  at  gradually  increasing  angles  and 
at  time-varying  crack  tip  speeds.  The  details  of  that  analysis  have 
been  completed  [19].  Publication  [19]  also  includes  a  study  of  the 
actual  conditions  for  crack  kinking  by  the  use  of  the  criterion  of  the 
balance  of  rates  of  energies. 
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